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Abstract 

The electric quadrupole and magnetic octupole moments of the light decuplet 
baryons are calculated in the framework of the light cone QCD sum rules. The 
obtained non-vanishing values for the electric quadrupole and magnetic octupole mo- 
ments of these baryons show nonspherical charge distribution. The sign of electric 
quadrupole moment is positive for fi^, S*" and negative for S*"^, which corre- 
spond to the prolate and oblate charge distributions, respectively. A comparison of 
the obtained results with the predictions of non-covariant quark model which shows 
a good consistency between two approaches is also presented. Comparison of the ob- 
tained results on the multipole moments of the decuplet baryons containing strange 
quark with those of A baryons shows a large SU(3) flavor symmetry breaking. 
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1 Introduction 



Detailed study of the electromagnetic properties of baryons, such as electromagnetic mul- 
tipole moments and electromagnetic form factors, can give essential information about the 
nonperturbative structure of QCD. These multipole moments arc related to the spatial 
charge and current distributions in baryons. Therefore, calculating these parameters could 
provide valuable insight on the internal structure as well as the geometric shape of baryons. 
The dominant elastic form factors of decuplet baryons are the charge Geq and magnetic 
dipole Gjvfi- The subdominant form factors are the electric quadrupole and magnetic 
octupolc GM:i (^11 these form factors are defined below). Note that, at = 0, the form 
factors Gmi, Ge2 and Gm-^ give the magnetic dipole /x^, electric quadrupole Qb, and the 
magnetic octupole moments Ob-, respectively [1]. The size of the higher multipole moments 
Qb and Ob provide information about the deformation of the baryon and its direction. 

Few words about the experimental prospects for measurement of the multipole moments 
are in order. There are two types of transitions for studying the multipole moments of 
the ground state decuplet baryons: diagonal transitions between them and off diagonal 
transitions between the decuplet and octet baryons i.e., A ^ A^, E* — E, E* ^ A and 
S* — > S. The couplings of diagonal decuplet-decuplet-photon transitions, obviously, can 
be measured only by virtual photon exchange. The magnetic moment of A"*" has been 
measured via 7p — > Tr^'j'p reaction [2]. However, measurement of the electric quadrupole 
by studying the diagonal transition is practically hopeless. This is due to the fact that the 
electric quadrupole operator is T-odd quantity and matrix element of this operator between 
the same initial and final states is equal to zero. Therefore, for the experimental study of 
the electric quadrupole moment, the suitable place is off-diagonal transitions. For example, 
the E2 transition can be measured in reaction octet baryons + X —>■ decuplet baryons + X 
[3], where X is heavy nucleons and also kaon photoproduction experiments 7p — > K+ 
decuplet K+ octet +7 [4] . Analysis of the electron-proton and photon-proton scattering 
experiments leads to a nonzero quadrupole moment of p — > A"*" transition [5]. 

There are large number of works in literature which are devoted to the investigation 
of the magnetic moment of hadrons, but unfortunately relatively little is known about the 
other multipole moments. Therefore, further detailed analysis is needed in studying higher 
multipole moments of the hadrons. Since obtaining direct experimental information about 
the electromagnetic multipole moments of these baryons is very limited, the theoretical 
studies play important role in this respect. The electric quadrupole and magnetic octupole 
moments of the A baryons have been calculated within the frame work of the light cone 
QCD sum rules (LCSR) in [6]. It is obtained that both quadrupole and magnetic octupole 
moments have nonzero values and negative sign, for example, for A+, implying that the 
quadrupole and octupole moment distributions of A"*" are oblate and have the same geomet- 
ric shape as the charge distribution. The same result has also been obtained by analyzing 
the quadrupole and octupole moments of A baryons in spectator quark model [7]. These 
multipole moments for A baryons have also been discussed in constituent quark model with 
configuration mixing but no exchange currents (impulse approximation), and constituent 
quark model with exchange currents but no configuration mixing [8]. 

Present work is devoted to the calculation of the electric quadrupole and magnetic 
octupole moments of the decuplet baryons in the framework of the light cone QCD sum 
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rules. As has already been noted, these multipole moments of the A baryons have been 
calculated in [6] in the same framework. Here, we extend the calculation of the multipole 
moments to the other members of the decuplet spin 3/2 baryons, i.e., S*"*"'*^'", and 
Q*~ . Note that, the magnetic octupole moments of these baryons have been calculated 
in non-covariant quark model (NCQM) [9]. Recently, the electromagnetic form factors of 
decuplet baryons have been calculated in lattice QCD in [10]. Here, also we stress that the 
magnetic moments of the decuplet baryons have been studied in [11] whithin light cone QCD 
sum rules. The main difference between the present study and [11] is that, in the present 
work we calculate additional form factors corresponding to different kinematical structures 
which are related to the higher multipole moments such as quadrupole and octupole. The 
outline of the paper is as follows: in section II, the light cone QCD sum rules for the 
electromagnetic form factors are obtained in LCSR. Section HI encompasses the numerical 
analysis of the form factors, a comparison of the results with the predictions of the other 
approaches and discussion. 

2 Light cone QCD sum rules for electric quadrupole 
and magnetic octupole moments of the Decuplet 
baryons 

For study the properties of hadrons in the sum rule formalism, the main working tool is 
the correlation function. To calculate the multipole form factors of the decuplet baryons, 
we consider the following correlation function: 

T^u = ij d'xe'P^O I T{r]^ix)UO)} I 0)^ , (1) 

where rj^ is the interpolating current for the decuplet baryons, and 7 denotes the elec- 
tromagnetic field. In the sum rule method, the above-mentioned correlation function is 
calculated in two different ways: on the phenomenological or physical side, it is saturated 
by a tower of baryons with the same quantum numbers as their interpolating current. On 
the QCD or theoretical side, it is calculated using the operator product expansion (OPE), 
where the short- and long-distance quark-gluon interactions are separated. The former is 
calculated using QCD perturbation theory, whereas the latter are parameterized in terms of 
the light-cone distribution amplitudes of the photon in light cone version of QCD sum rules. 
The electromagnetic form factors are determined by matching these two representations of 
the correlation function. 

First, let us calculate the physical part of the correlation function. By isolating the 
contributions of the ground state baryons from Eq. ([T]), we obtain 

^ I + . . ■ , (2) 

P2 — PI — 

where pi = p + q, P2 = P and q is the momentum of a photon. The dots mean contributions 
of the higher states and continuum. 



2 



It follows from Eq. ([2]) that, for calculation of the physical part, we need to know the 
matrix element of the interpolating current between the vacuum and the decuplet baryon 
state as well as transition matrix element, (-8(^2) | B{pi))^. The (0 | ?7/x(0) | B{p,s)) is 
defined in terms of the residue of the corresponding decuplet baryons, as: 



{0\r],iO)\Bip,s)) = \BU^{p,s), 



(3) 



where u^{p^s) is the Rarita-Schwinger spinor. The transition matrix element {B[p2) 
B(jpi))^ can be parameterized in terms of four form factors as [7, 12, 13]: 
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+ 



(2m 



LL V 



B 



2mB 

Uu{pi), 



F2g^'' + F, 



{2mBy 



^4 



(4) 



where e is the polarization vector of the photon and Fi are form factors as functions of 
transfer momentum square q^ = (^1—^2)^- For obtaining the expression for the correlation 
function from physical side, summation over spins of the spin 3/2 particles is performed 
using 



^Uf,{p,s)u,y{p, s) = i^ + mB) 



1 



'^Pt^Pu Pi^lu - Pull, 



3m| 



3mB 



(5) 



In principle, using the above equations, we can obtain the final expression of the physical 
side of the correlation function, but we come across with two difficulties, a) Not only spin 
3/2, but spin 1/2 particles contribute to the correlation function, i.e., the matrix element 
of the current 77^ of the spin 3/2 particles between vacuum and spin 1/2 states is nonzero. 
This matrix element in general form can be written as 



(0 I r/^(0) I Bip, s = 1/2)) = (Ap^ + B^,)u{p, s = 1/2). 



(6) 



Using the condition 7^77'^ = 0, one can immediately obtain that B = —jm. b) All Lorentz 
structures are not independent (for more details, see [14]). 

In order to eliminate the unwanted spin 1/2 contributions and obtain only independent 
structures, the ordering procedure of Dirac matrices are applied and in the present work, 
we choose it as 7^ ^ ^ ^j^. After this ordering procedure, we obtain the final expression of 
the physical side of the correlation function as follows: 



T, 



A 



1 



2(e ■ p)gi,u i>Fx + — (£ ■ p)g^y i> iF2 
rriB 



+ 



^ {pI - m\){pl - m\) 

' ^ ■ p)%(lv iiF-i + - — 2"(^ ■ P)%'iv iF/^ + other independent structures 



2m| 



4m| 



+ structures with 7^ at the beginning and '^y at the end or which are 
proportional to p^p, or p^y 



(7) 
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For calculation of the four form factors, wc need four structures. We will choose the 
structures {e.p)g^j^y {e.p)g^i, t/) ^, {b-p)^^^^ i/> and [E.p)q^qy ^ for determination of the form 
factors Fi, F2, F3 and F4, respectively. In the experiments the multipole form factors, Ge^ 
(charge), Gmi (magnetic dipole), (electric quadrupole) and Gm^ (magnetic octupole) 
are usually measured. Therefore, we need relations between two sets of form factors. The 
multipole form factors are defined in terms of the form factors Fi{q^) as [7, 12, 13, 15]: 

GeM') = [m') - xF2{q')\ {l + - [F,{q') - xF,{q^)\ | (1 + x) , 

GmM') = [H<l') + Hq')] (1 + [^3(g') + F,{q')] x{l + x) , 

GeAq') = [Fi{q') - xF^iq')] - \ [F,{q') - xF,{q^)\ (1 + x) , 

GmA^) = [F,{q') + F,{q')\ - \ {F^{q^) + FM^)\ (1 + x) , (8) 

where x — —q^/im\. At q^ = 0, we obtain 

G'mi(O) - Fi(0) + F2(0) , 
GeM = Fi(0) - 1^3(0) , 

Gms{0) = Fi(0) + F2(0)-^[F3(0) + F4(0)] . (9) 

The magnetic dipole /ib, the electric quadrupole Qb, and the magnetic octupole Ob mo- 
ments are defined in terms of these form factors at = in the following way: 

g 

f^B = 7^ Gmi(0) , 

zmB 

g 

Qb — — 2-Ge2{0) ■, 

Ob = i^GmM- (10) 

The QCD side of the correlation function, on the other hand, can be calculated by the 
help of the OPE in deep Euclidean region where -C and (p + g)^ -C 0. For this aim 
we need to know the explicit expressions of the interpolating currents of the corresponding 
baryons. The interpolating currents for decuplet baryons are [16] 



y*0 






















= 71"' 



y|e"^^[(M'^^C7^d'')s^ + {d'^^G'y^s'X + (5'^^C7^«'')d1 , 



Q*— 1 y*+ / X 

(11) 

where a, b and c are color indices and C is the charge conjugation operator. After contracting 
out all quark pairs in Eq. ([T]) using the Wick's theorem, we obtain the following expression 
for the correlation function of the — > S*°7 transition in terms of the quark propagators: 

I Qc'*'- Q"^<^' ^ ofec' 1 Qca' Qfbb' oac' , Qcb' Qiaa' nbc' 

+ rr(7,^f 7.^^Vf } I 0), (12) 

where S' = CS^C and Su,d,s are the hght quark propagators. The correlation functions 
for other transitions can be obtained by the replacements mentioned in Eq. (fTTj) . The 
expression of the light quark propagator in the external field is calculated in [17, 18]: 

- igs I du l -i^'t ^ G^^{ux)a^'' - ux^G^yiux)^ 
Jo I ^^'^ ^ 



-^^G,,{ux)a>^^ 



327r 



In ( + 27^ 



(13) 



where A is the scale parameter, and following [19], we choose it at the factorization scale 
A = 0.5 GeV — 1.0 GeV . The correlation function contain three pieces: a) Short distance 
contributions, b) "Mixed" contributions, c) Large distance contributions when a photon is 
radiated at long distance. 

Different terms in Eq. f|T3|) give contributions to the different pieces of the correlation 
function. The short distance contributions can easily be obtained from Eq. f[T^ by replacing 
one of the propagators by. 



S'-^p = y d'ySf^^^x - y) 4Sf''^\y) \ , (14) 
where S^'^'^'^ is the light quark propagator given as. 



= ^ - ■ (15) 



i ^ rriq 
2t:'^x^ An^x'^ 

and two other quark propagators are replaced by the free quark propagator. 

In the "mixed" contributions photon interacts with quark fields perturbatively. 

Therefore, one of the quark propagators is replaced by Eq. dHj), two other propagators 
either both are replaced by 

S^'f, = -\rr,q\T,U , (16) 
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which both can form quark condensates, or one of them is replaced by Eq. (fT6|) and second 
one by the free quark propagator. In Eq. (flGl) . Tj is the full set of Dirac matrices. 

The large distance contributions can be obtained from Eq. ( fT2l) by following replace- 
ments: One of the quark propagators is replaced by Eq. ( fT6l) and a photon interacts 
with the quark fields at large distance, i.e., the matrix elements of the nonlocal operators 
q{xi)Tq'{x2) and q{xi)G^^Tq'{x2) appear between the vacuum and the vector meson states, 
which is parameterized in terms of photon distribution amplitudes (DA's). Two other prop- 
agators are both replaced by free quark propagator, or one of them is replaced by free quark 
propagator and second one is replaced by Eq. (fT6l) and then it interact with QCD vacuum, 
i.e., it forms a quark condensate, or both of propagators are replaced by Eq. ( fT6i) and then 
they form quark condensates. 

Using the expressions of the light propagators and the photon DA's and separating the 
coefficient of the structures mentioned before and applying double Borel transformation 
with respect to the variables pi = and pi = {p + q)'^ to suppress the contributions of 
the higher states and continuum, sum rules for the form factors Fi, F2, F3 and F4 are 
obtained. The explicit expressions of the sum rules for these form factors are given in 
the Appendix-A. From the expressions of the form factors it is clear that, to obtain form 
factors, we need to know the explicit expressions of residues of the corresponding baryons. 
The explicit expressions for these residues are given in [20,23]. 

3 Numerical analysis 

Present section is devoted to the numerical analysis for the, electric quadrupole and mag- 
netic octupole moments of the light spin 3/2 baryons. The values for input parameters used 
in the analysis of the sum rules for the Fi, F2, F^ and F4 are : (^^(l GeV)) = {dd{l GeV)) = 
-(1.65 ± 0.15) X 10-2 GeV^ [21], {ss{l GeV)) = 0.8{uu{l GeV)), m,(2 GeV) = (111 ± 
6) MeV at Aqcd = 330 MeV [22], ml{l GeV) = (0.8 ± 0.2) GeV^ [23] and /g^ = 
—0.0039 GeV^ [24]. The value of the magnetic susceptibility is taken to be x(l GeV) = 
—3.15 ± 0.3 GeV^^ [24]. As has already be noted, the main input parameters in light cone 
sum rules are the DA's. The explicit expression of the photon DA's are given in [24]. 

The sum rules for the electromagnetic form factors also contain two auxiliary parame- 
ters: Borel mass parameter and continuum threshold Sq. The physical quantities should 
be independent of these parameters. Therefore, we look for a region for these parameters 
such that the electromagnetic form factors are independent of them. The working region for 
are found requiring that not only the contributions of the higher states and continuum 
should be less than the ground state contribution, but the highest power of l/M^ be less 
than say 30*^/0 of the highest power of M^. These conditions are satisfied in the regions 
1.1 GeV^ < M2 < 1.6 GeV^, 1.2 GeV^ < < 1.7 GeV^ and 1.4 GeV^ < < 2.4 GeV^ 
for S*, S* and Q* baryons, respectively. In the numerical analysis, Sq = {ms + 0.5)^ GeV^ 
has been used for value of the continuum threshold. 

Our final results on the electric quadrupole Qb and magnetic octupole Ob moments of 
decuplet baryons are presented in Table 1. The quoted errors in Table 1 can be attributed 
to the uncertainties in the variation of the Borel parameter M^, the continuum threshold 
So, as well as the uncertainties in the determination of the other input parameters entering 
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the sum rules. A comparison of our predictions on magnetic octupole moment with the 
results obtained in NCQM is also presented in Table 1. The results for magnetic octupole 





Quadrupole Q^ffTi^) 


Octupole 


O(fm^) 




Present Work 


Present Work 


NCQM [9] 


n- 


0.12±0.04 


0.016±0.004 


0.003-=-0.012 




0.03±0.01 


0.013±0.004 


0.008^0.012 




0.0012±0.0004 


-0.001±0.0003 


0.000^0.002 


E*+ 


-0.028±0.009 


-0.015±0.005 


-0.004^-0.012 




0.045±0.015 


0.020±0.006 


0.005^0.012 


"*0 


0.0025±0.0008 


-0.0014±0.0005 


0.000^0.002 



Table 1: Results of the electric quadrupole moment (in units of jvr?) and magnetic octupole 
moments (in units of /m^) of the decuplet baryons. 

moments show a good consistency between our predictions and those of the NCQM [9]. As 
has already been noted, the electromagnetic form factors of the decuplet baryons have been 
calculated at 7^ in [10], so we can not compare our results with theirs. However, the 
order of magnitude of our results are in good agreement with their predictions at low 
Comparison between our results on electric quadrupole and magnetic octupole moments 
of E*+'~, S*~, Or and the predictions of [6] for A baryons, shows a large SU(3) flavor 
symmetry braking. This violation is larger for S*~ and Or baryons which contain two and 
three strange quarks, respectively. In the case of the strange baryons, the results are very 
sensitive to the strange quark mass. This sensitivity together with the different working 
regions of Borel mass parameter, , and continuum threshold, sq lead to the different 
values of multipole moments for E*, 5*, O baryons. 

In conclusion, the electric quadrupole and magnetic octupole moments of decuplet 
baryons were calculated in the framework of the light cone QCD sum rules. We obtained 
non- vanishing values for the electric quadrupole and magnetic octupole moments of these 
baryons which mean nonspherical charge distribution. The sign of electric quadrupole mo- 
ment is positive for ^2", S*^ and negative for E*"*", which correspond to the prolate 
and oblate charge distributions, respectively. The obtained results are in good consistency 
with the predictions of the non-covariant quark model. Comparison of the obtained results 
on the multipole moments of the decuplet baryons containing strange quark with those of 
A baryons presents a large SU(3) flavor symmetry breaking. 
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Appendix A 



In this appendix, we present the sum rules for the form factors, -Fi(O), -F2(0), -^3(0) and 

^4(0). 



Fi{q^ = 0) 



1 



407r4 



(e„ + ed + e,)M^ 



67r2 



2(e<i + es){uu) + 2(e„ + es)((ici) - (e^ + ed){ss) 



1 



+ 



+ 



367r2M2 
1 

547r2M2 
1 

1447r2M4 

1 2 
-mnrris 



A2 

ms{giG^){eu{dd) + ed{uu)) ( 7^ + In — 



nisig^G ){eu{dd) + ed{uu)) 



9M2 



ml{eu{dd) {ss) + ed{uu){ss) + es{uu) {dd)) 



mQirisig^G ){e^{dd) + ed{uu)) 



(9ed + 10e,)(MM) + (9e„ + 10e,)(rf(i) - 4(e„ + ed){ss) 



+ -{eu{dd){ss) + ed{uu){ss) + es{uu){dd)) | , 

F,{q' = 0) = !^e-^-/M^ / ^—-:m, ( + In 

X {[3M^ + 27r^Ma{uoWsG^){eu + e^) + 24e,M6} 



1 



+ 



1 

144^ 



M'^[ms{3eu + Se^ + lle^) + 87r^x(e„(MM) + ed((i(i) + es{ss))(l)^{uo)] 
M2{6[(e„ + ed)(ss) + (e„ + e,)(dd) + (e<i + es){uu)] 



+ {e^{uu) + ed{dd) + e,(ss)) (^3A(mo) - 4[i2(5, 1) + ^2(5, 3 - 4.v) 
+ i2(T2, 1 - 2^;) + 2i2{Tz, l-2v)- i2{%, 1 - 2v) + Sisih^)]^ 
- 3ms/3^(e„ + erf)^„(Mo)} 

ms{ss){eu{uu) + ed{dd))[i2{S, 1) + ^2(5, 3 - 4^) + i2(^, 1 - 2v) 



54M2 



+ 2i2(r3, l-2v)- 1 - 2t;)] + 



2{eu{dd) + ed{uu)){ss) 



36M2"' 

+ {eu{uu) + ed{dd)){ss)A{uo) + 16(e„ + ed){uu){dd)i2,{h^) 

1 „ , o„o 1 



+ 



+ 



msfsyig^G )(e„ + ed)V'a(^^o) + 



(A.1) 



ll/37[(ed + es)(w'w) + (e„ + es)((ici) + (e„ + ed){ss)]ipa{uo)^ 



mQ< 2Ams{ss)x{eu{uu) + ed{dd))(f)^{uo) 



morris 



[eu + ed){uu){dd) - 2{eu{uu) + ed((i(i))(ss) U3(/i^) 
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-ms{glG'^){eu + e^) - 7^7r^"^o[(ed + es)('UM) + (e„ + e5)((i(i) 



11527r4 /V " ■ 432^2 



+ (e„ + ed){ss)] + ^{ - 2msx{eu{uu) + ed{dd)){ss)4>^{uQ) + fsjUed + es){uu) 

+ (e„ + e,)((id) + (e„ + ed)(ss)]V'aK) } (A.2) 

- ^MV37(e„ + ed + e,) ([2i2(A 5 - 4v) + 4i2(V, 1 - 2^;) - V'aK)]) 



+ 3gJ^4 m,(e^(m) + ed((i(i)) (4M^[5zi(Ti + T2, 1) - 3zi(T3 + T4, 1)] 



A2\ 
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1 4 - 

+ 54M%2 ^°^'^^^^^^" ^ ^'^^ ^ 27M2 f^^" ^ ed){uu){dd) + (e„ + e,) (iZii) (ss) 

+ (ed + e,)((io?)(ss)][?i(3Ti + 4T2 - T4, 1) + Qi^{h^)] 

- ^^msf3j{ss){eu + ed)[i2iA, 5 - 4^;) + 2i2(V, 1 - 2^;) - 3V'a(tio)] 

+ 215^4^2 + ed{dd))[zi{3T, + 4r2 - %, 1) + 4T3(/i-,)] 

2 / - _ ~ 

- ^Y]f^"^o(^0[(e« + ed){uu){dd) + (e„ + es){uu){ss) + (e^ + es)(dd)(ss)]i3(/i^) 

+ rrisfs^ieu + ed){ss)tlJa{uo) 
1 



3(e„ + ed)(ss) - 8{eu{uu) + ed{dd))[h{2Ti + ?,% - 2%, 1) 



- sTsl/i^)]) (A.3) 



- Yii;;^^'-^^^^''" + '''^ + ""^^ (4[4^2(A 1 + + 4^2(V, 1 - v) +?3(^.)] - 3V^,(mo) 



4 _ - 

+ l^^f^ [(^« + ^'i) i^d) + (e„ + e^) {uu) (ss) 
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- ^Ym^"^°(^^^*^^" ed){uu){dd) + (e„ + es){uu){ss) + (e^ + es)((irf)(ss)]i3(/i^) 
+ rrisf 3j{eu + ed)(ss)V'a(iio)) 

- T^TTT'si^eu + ed)(ss) - 4(e„(MM) + ed{dd))[4ii{Ti - T3, 1) 



where, the Borel parameter is defined as = MfM^/M^ + M| and = Mf/{Mf + 
M|). Since the masses of the initial and final baryons are the same, we have set Mf = M| 
and Uq = 1/2. The continuum subtractions have been made via M^" — >■ M'^^En{x), where 
E^ix) = 1 - e-- E:=o' f witji X = s,/M\ 

The functions ^3 and is are also defined as 



108M2 
1 



ms{glG''){e^{uu) + eSdWk{T2-%,v)+H%)] 



+ 



108M%2 
1 





(A.5) 



where k — aq-\- agV. 
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